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A complete theoretical study is presented for the disappearance of flow, for the first time, by 
analyzing 15 reactions with masses between 47 and 476 units. We demonstrate that the effect of 
nucleon-nucleon cross-section reduces to insignificant level for heavier colliding nuclei in agreement 
with previous studies. A stiff equation of state with nucleon-nucleon cross-sections cr=35-40 mb is 
able to explain all the measured balance energies within few percent. A power law {onA'^) is also 
given for the mass dependence of the disappearance of flow which is in excellent agreement with 
experimental data. 
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I. INTRODUCTION 

The heavy-ion collisions at intermediate energies pro- 
vide a rich physical insight into the reaction dynamics. 
One has measured (and/or predicted) several new phe- 
nomena that may shed light on the nature of hot and 
dense nuclear matter formed during a collision. In ad- 
dition, one also hopes to understand the nature of nu- 
clear interactions in medium. The prediction of collec- 
tive transverse flow by the hydrodynamical model was 
a very important step towards the understanding of ex- 
cited nuclear matter 0. The collective transverse flow 
was found to be very sensitive towards different signals of 
excited nuclear matter. Apart from the transverse flow, 
one has also proposed e.g. differential flow .2|, and el- 
liptic flow jsj- All these quantities are assumed to be 
sensitive towards the (nuclear matter) equation of state 
and/or nucleon-nucleon (nn) cross-section=> the ultimate 
goal of the intermediate energy heavy-ion collisions. One 
should however, keep in the mind that the reaction dy- 
namics depends also on the incident energy as well as 
on the impact parameter of the reaction. At low inci- 
dent energies, the dynamics is governed by the attractive 
mean field whereas the repulsive interactions decide the 
fate of a reaction at higher incident energies. Naturally, 
the effect of nn collisions decreases with decrease in the 
incident energy. The dominance of the attractive mean 
field (at low incident energies) may prompt the emission 
of particles into backward hemisphere whereas if nn scat- 
terings dominate, the particle emission is likely to be in 
the forward hemisphere. Therefore, while going from the 
low incident energy to higher energy, the attractive in- 
teractions may be balanced by the repulsive interactions, 
resulting in the net zero flow (i.e. the disappearance of 
flow) . The energy at which the flow disappears is termed 
as the balance energy P| . 

During the last few years, extensive efforts have been 
made to measure and understand the disappearance of 

flow iEiiiiiiiainiiiiiiiiiiiiiiiiiia. 

One has measured the balance energy (Ehni) in ^^C-|-^^C 
H, 20Ne+2Ul_i 36Ar+27Al Jl7|, 40Ar+27AlJ|, 
40Ar+45Sc 11 a d, 40Ar+5iV [J^ 642^+27 ai Fq^ 
40Ar+58Nj 64Zn+48Ti [l^, ^^Ni-h^SNi [Hill El, 



58pg_^58pg ^1 642n+58]s^i [T3L j^r-^^^^b [1 111, 

93Nb+93Nb i29xe+"8snlll, i39La+i39La [3, 

and i^Uu-h^^^Au EHIlIll systems. The very recent 
and accurate measurement of the balance energy Ebai in 
197Au+197Au [ll m has generated a renewed interest 
in the field • Interestingly, most of the reported reac- 
tions are symmetric in nature. It should be kept in the 
mind that the reaction dynamics depends also upon the 
asymmetry of the reaction |T^ . All the above mentioned 
measurements were for the central collisions only. Some 
measurements 0,11,0,0,0,01 however, also took the 
impact parameter dependence into account. As noted in 
ref. [ill, the E^ai for heavier nuclei shows a little depen- 
dence on the impact parameter whereas a large variation 
in the Ebai can be seen for lighter colliding nuclei |^ ■ 
The possible cause could be the fact that the disappear- 
ance of flow for heavier nuclei occurs at a much lower 
incident energy compared to lighter nuclei (e.g. the mea- 
sured Ebai for i^'^Au-h^^^Au is 42±3±1 MeV/nucleon 
whereas it is lllilO MeV/nucleon for ^°Ne+^^Al p). 
In the (nearly) absence of nn collisions at low incident 
energies, a little dependence should occur on impact pa- 
rameter |l9l| . Some attempts also exist in the literature 
where enhancement in the E^ai with neutron content was 
found experimentally/theoretically |T^l20ll2ll|. 

The above findings reveal the measurements of balance 
energy in more than 15 systems ranging from very light to 
heavy nuclei. As a result, a power law behaviour (ocA^) 
has also been reported for ,^5, 12, 16, 17]. Earlier the 
power law parameter "r" was supposed to be close to - 
1/3 (resulting from the interplay between the attractive 
mean field and repulsive nn collisions) @ whereas the 
recent measurements suggest a deviation from the above 
mentioned power law (l2l lla | . 

Various theoretical attempts have been made to under- 
stand the vanishing of nuclear flow. Most of these are, 
however jUsing the Boltzmann-Uehling-Uhlcnbeck model 
II S i i & iO, ,12, il, U, 16, 21, 22, 23, 24, M- 
Some attempts are also reported in the literature where 

Suantum Molecular Dynamics (QMD) model was used 
i 20, 27, 28, 29, 30, 31). Different theoretical attempts 
considered either a stiff or soft equation of state along 
with variety of nn cross-sections. Interestingly, out of 
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all these studies, only a couple of attempts exist where 
mass dependence of the disappearance of flow was dis- 
cussed la [p, E ^lM, Ha E3- a careful analysis 
shows that in refs. ly, '26'| a total mass <200 was consid- 
ered for the power law studies whereas in refs. 30, .sj 
only heavier masses >175 were analyzed. There the Ebai 
for heavier nuclei was found to scale as approximately 
30, 31] whereas lirfiter and medium mass nuclei 
follow A-i/3 dependence 13 [23. Refs. [Ulil included 
for the first time, a larger mass range 63/47<A<394. 
However, even in these studies, only six systems were 
taken into account (note that the power law fit was 
made for those systems which were recorded in the 47r 
NSCL experiment. The data from the GANIL exper- 
iments, however, was not taken into account for fit- 
ting). The earlier mass dependence by the same group 
13 also considered five systems only. We here present, 
for the first time, a complete study of the mass de- 
pendence in the disappearance of flow by analyzing as 
many as fifteen central collisions ranging from ^°Ne-|-^^Al 
to i^'^Au-l-^'^'^Au/^^^U-l-^^^U. The corresponding experi- 
mental measurements were re p orte d by several different 

groups iillSi|iaiii|il|iIlil|il|illi3. 

This study, therefore, will present a unified analysis 
of the balance energy irrespective of the experimental 
setup/group. The present study is carried out within the 
framework of the Quantum Molecular Dynamics (QMD) 
model [lllliliSlHllillillilslllilsllll. The de- 
tails of the model are given in section 2. The results and 
discussion is presented in section 3 and we summarize 
our results in section 4. 
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where the total time derivative includes the derivatives 
with respect to the parameters. The time evolution is 
obtained by the requirement that the action is stationary 
under the allowed variation of the wave function 
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*]dt = 0. 
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If the true solution of the Schrodinger equation 
is contained in the restricted set of wave function 
(j)a {xi, Xa,Pa) j this variation of the action gives the ex- 
act solution of the Schrodinger equation. If the parame- 
ter space is too restricted, we obtain that wave function 
in the restricted parameter space which comes close to 
the solution of the Schrodinger equation. Performing the 
variation with the test wave function (2), we obtain for 
each parameter A an Euler-Lagrange equation; 



(6) 



For each coherent state and a Hamiltonian of the form. 



the Lagrangian and the 



Euler-Lagrange function can be easily calculated |33 
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II. THE MODEL 
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We describe the time evolution of a heavy-ion reaction 
within the framework of Quantum Molecular Dynamics 

(QMD) model ^^Qg,2L2^,22^M,M^^M,M 
which is based on a molecular dynamics picture. Here 
each nucleon is represented by a coherent state of the 
form 
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(j)a{xi,t) = ( — 



(1) 



Thus, the wave function has two time dependent param- 
eters Xa and Pa- The total n-body wave function is as- 
sumed to be a direct product of coherent states: 



4> = 4>aixi,Xa,Pa,t)<j)f3ix2,Xi3,Pi3,t). 



(2) 



where, antisymmetrization is neglected. The initial val- 
ues of the parameters are chosen in a way that the ensem- 
ble {Ax+Ap) nucleons give a proper density distribution 
as well as a proper momentum distribution of the projec- 
tile and target nuclei. The time evolution of the system 
is calculated using the generalized variational principle. 
We start out from the action 
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Thus, the variational approach has reduced the n-body 
Schrondinger equation to a set of Gn-different equations 
for the parameters which can be solved numerically. If 
one inspects the formalism carefully, one finds that the 
interaction potential which is actually the Bruckner G- 
matrix can be divided into two parts: (i) a real part and 
(ii) an imaginary part. The real part of the potential 
acts like a potential whereas the imaginary part is pro- 
portional to the cross-section. 

In the present model, interaction potential comprises 
of the following terms: 



(10) 



Vioc is the Skyrme force whereas Vcoui and Vvuk define, 
respectively, the Coulomb and Yukawa terms. The ex- 
pectation value of these potentials is calculated as 

xd^r^d^rpd^Pad^pp, (11) 
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Where (p^ , , t) is the Wigner density which corre- 
sponds to the wave functions (eq. 2). If we deal with the 
local Skyrme force only, we get 
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Here A, B and C are the Skyrme parameters which 
are defined according to the ground state properties of a 
nucleus. Different values of C lead to different equations 
of state. A larger value of C (=380 MeV) is often dubbed 
as stiff equation of state. 

A number of attempts exist in the literature which 
study the nature of e quat ion of state. Following refs. 

i II i i i 113 MM El m m m na lii, we 

shall also employ a stiff equation of state through out the 
present analysis. It should also be noted that the success 
rate is nearly the same for stiff and soft equations of state. 
Further, it has been shown in rcfs. uSi, .22, ,24. ^26] that 
the difference between Ebai using a stiff and soft equation 
of state is insignificant for central heavy-ion collisions. 

The imaginary part of the potential i.e. the nn cross- 
section has been a point of controversy. A large number 
of calculations exist where a constant and isotropic cross- 
section is used. Following refs. H [H IH [H Hi H [H 
I35I l3^ , we also use constant energy independent cross- 
section. As shown by Li , most of the collisions below 
100 MeV/nucleon happen with nn cross-section of 55 mb 
strength. Keeping the present energy domain into mind, 
the choice of a constant cross-section is justified. It has 
also been also shown by Zheng et al. 3] that a stiff 
equation of state with free nn cross-section and a soft 
equation of state with reduced cross-section yield nearly 
the same results. For comparison, we shall also use an 
energy dependent cross-section as fitted by Cugnon [s^l 
(labeled as Cug) as well as a medium dependent cross- 
section derived from G-matrix js^] (denoted by GMC). 



III. RESULTS AND DISCUSSION 

Using a stiff equation of state along with different nn 
cross-sections, we simulated the above mentioned reac- 
tions for 1000-3000 events in each case. The reactions 
were followed till transverse flow saturates. The satu- 
ration time varies between 150 fm/c (for lighter collid- 
ing nuclei such as ^"Ne-l-^^Al) and 300 fm/c (for heav- 
ier colliding nuclei such as ^^^Au-|-^^^Au). In particu- 
lar, we simulated '^°Ne+^'^A\ at b=2.6103 fm, ^^At+^'^A\ 
at b=2 fm, 40Ar-h27Al at b=1.6 fm, ^OAr+^SSc at 
b=3.187 fm, "OAr-h^iy at b=2.442 fm, ^^Zn+^'^A\ at 
b=2.5 fm, ^OAr-h^^Ni at b=0-3 fm, ^^Zn+^Sxi at b=2 



fm, 58Ni-H58Ni at b=2.48 fm, ^'^Zn+^^Ni at b=2 fm, 
ssKr-h^^Nb at b=4.07 fm, ^^Nb-f^^Nb at b=3.104 fm, 
i29Xe-h^"*Sn at b=G-3 fm, "^La+^SLa at b=3.549 fm 
and ^9^Au-|-"'^^^Au at b=2.5 fm. The choice of impact 
parameter is based on the experimentally extracted in- 
formation [5-17]. The above reactions were simulated 
at incident energies between 30 MeV/nucleon and 150 
MeV/nucleon depending upon the mass of the system. 
Naturally, a lower energy range was used for heavy nu- 
clei whereas a higher beam energy was needed for lighter 
cases. The reactions were simulated at different fixed in- 
cident energies and a straight line interpolation was used 
to find the balance energy Etai- 

There are several methods used in the literature to de- 
fine the nuclear transverse flow. In most of the studies, 
balance energy is extracted from (jj^/A) plots where one 
plots (px/A) as a function of Ycm/Yheam- Using a lin- 
ear fit to the slope, one can define the so-called reduced 
flow (F). Naturally, the energy at which the reduced flow 
passes through zero, is called the balance energy. Al- 
ternatively, one can also use a more integrated quantity 
"directed transverse momentum (p!^*'')" which is defined 
as mililMillsllslIll: 



(14) 



dir \ 



A 



where Y{i) and Pxii) are the rapidity distribution and 
transverse momentum of the ith particle. In this defini- 
tion, all rapidity bins are taken into account. It, there- 
fore, presents an easier way of measuring the in-plane 
flow rather than complicated functions such as the (px/A) 
plots. It is worth mentioning that the balance energy is 
independent of the nature of emitted particle Fur- 
ther, the apparatus correction and acceptance does not 
lay any role in calculating the energy of vanishing flow 



In fig. 1, we display the flnal state transverse momen- 
tum (px/A) as a function of the rapidity which is defined 
as: 



^' 2 E(i)-p,(f)' 



(15) 



where E(i) and Pz(j) are, respectively, the total energy 
and longitudinal momentum of ith particle. The upper 
parts are for ^°]<ie+^'^Al (at 150 fm/c) and "OAr-f^SSc (at 
150 fm/c), whereas the bottom parts are for ^^^La-l-^^^La 
(at 300 fm/c) and i^^Au-h^^^Au (at 300 fm/c). The mid- 
dle part is for ^^Zn+^^ Ni (at 300 fm/c) and ^^Nb-l-^^Nb 
(at 300 fm/c). In all cases, the slope is negative at 
lower incident energies which changes to positive value 
at higher incident energies. Between these limits, the 
slope becomes almost zero at a particular energy. This 
zero slope energy is termed as balance energy. One also 
notices that a higher value of incident energy is needed 
in lighter cases to balance the attractive and repulsive 
forces. This energy decreases with increase in the mass 
of the system. 
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A look at fig. 2, where (p^*'") (instead of {px/A)) is 
plotted, depicts quite similar trends. Here (p^"") is dis- 
played as a function of the reaction time. The (p**^) dur- 
ing initial stage is always negative irrespective of the inci- 
dent energy. This shows that the interactions among nu- 
clei are attractive during the initial phase of the reaction. 
These interactions remain either attractive throughout 
the time evolution, or may turn repulsive depending on 
the incident energy. The transverse flow in lighter collid- 
ing nuclei saturates earlier compared to heavy colliding 
nuclei. One also sees a sharp transition from negative to 
positive flow in lighter nuclei. This transition is gradual 
when one analyzes the heavier nuclei. If one compares 
figs. 1 and 2, one finds that the disappearance of flow 
(where flow is zero) occurs at same incident energies in 
both cases showing the equivalence between {pxl A) and 
(pIJ"") as far as balance energy is concerned. The latter 
quantity is more useful since it is stable than former one. 
These findings are in agreement with refs. 28, 29]. 

It has been advocated by several authors that the study 
of disappearance of flow can shed l ight on the magnitude 
of nn cross-section |lllll[l3[ll|ll[ll|2l Hll IH 

llEll^ll^l^- To check this point with reference to mass 
dependence, we show in fig. 3, the final stage (p^"") as 
a function of the incident energy. The left panel is for 
20Ne-H27Ai^ 642n+58Nj aj^j i39lj^_^139lj^^ respectively, 

for top, middle and bottom whereas top, middle and bot- 
tom of right hand side displays the results, respectively, 
of 40Ar-f45Sc, 93Nb-H93Nb and ^^^Au-f ^^Uu. The exper- 
imental data are displayed by stars whereas the present 
results with cr=55, 40 and 20 mb are shown, respectively, 
by solid circles, open squares and solid triangles. The 
(p^*'') obtained with energy dependent cross-section due 
to Cugnon (labeled as Cug) and the one that takes the 
medium into account (i.e. the G-matrix) are marked by 
solid diamonds and inverted triangles, respectively. First 
of all, we see that the medium effects in nn cross-section 
do not play any role at these low incident energies. The 
results obtained with Cugnon energy dependent cross- 
section and G-matrix medium dependent cross-section 
are roughly the same for heavier colliding nuclei where re- 
actions are simulated at low incident energies. Some vis- 
ible differences, however, can be seen in the case of light 
colliding nuclei where incident energy is relatively high. 
One also sees a linear enhancement in the nuclear flow 
with increase in the incident energy. Further, the role 
of different cross-sections is consistent through out the 
present mass range. The largest cross-section gives more 
positive flow which is followed by the second larger cross- 
section. Interestingly, these effects depend on the mass 
of the system. If one looks at the reaction of ^"Ne-j-^^Al, 
one sees that the Ebai increases from 89 MeV/nucleon to 
244 MeV/nucleon when nn cross-sections is reduced from 
55 mb to 20 mb. Whereas the range of E^ai for the same 
cross-sections was 47-83 MeV/nucleon for ^'^Nb-l-^^Nb re- 
action. If one goes to still heavier nuclei, ^^^Au-|-^^^Au, 
the range of E^ai narrows down to 38-59 MeV/nucleon. 
In other words, a reduction in the cross-section by 64% 



percent yields a change of 155 MeV/nucleon in the case of 
^"Ne-l-^^Al reaction, whereas it is only 21 MeV/nucleon 
for the case of ^^^Au-|-^^^Au. Similarly looking at the 
curves of (t=55 and 40 mb, a reduction in the cross- 
section by 27% yields a difference of 30 MeV/nucleon in 
the Ehai for ^°Ne-|-^^Al reaction whereas nearly 4 MeV 
difference exists for the case of ^^^Au-|-^^^Au reaction. 
This result, which is in agreement with the findings of 
refs. 0, depicts that for heavier colliding nuclei, 
the Ebai is independent of the cross-section one is choos- 
ing. Further, the standard energy dependent nn cross- 
section (Cug) fails to reproduce the observed balance en- 
ergy in almost all the cases. However, a constant cross- 
section of 40 mb strength seems to be closer to the ex- 
perimental observed balance energy. This conclusion is 
supported by several other groups where a cross-section 
of 30-40 mb was used to reproduce the experimental data 
[lllllllillllllllllllalllllllll. This will be dis- 
cussed in detail in the following paragraphs. 

It has been argued in refs. |2^ [s^l that the flow at 
any point during the reaction can be divided into the 
parts emerging from the (attractive) mean field potential 
and (repulsive) nn collisions. Following |0, we decom- 
posed the transverse momentum into the contributions 
created by the mean field and two- body nn collisions. 
This extraction, which is made from the simulations of 
the QMD model, is done as following |28|]: Here at each 
time step during the collision, momentum transformed by 
the mean field and two- body collision is analyzed sepa- 
rately. Naturally, we get two values at each step which 
can be followed throughout the reaction. The total trans- 
verse momentum can be obtained by adding both these 
contributions. 

In fig. 4, we display the final state (p^"") decomposed 
into two parts i.e. into the mean field and two- body 
collision parts as a function of the incident energies for 
different colliding systems as reported in fig. 3. Again a 
linear enhancement in the flow with energy can be seen. 
Further, the contribution of the mean field remains at- 
tractive throughout the energy range whereas collision 
contribution is always repulsive. The balancing of both 
these contributions results in net zero flow. One should, 
however, keep in the mind that the contribution of the 
mean field potential may even turn repulsive at higher 
incident energies [37| . 

In fig. 5, we display the energy of vanishing flow (Ebai ) 
as a function of the combined mass of the system for fif- 
teen different colliding pairs that range from ^°Ne-|-^^Al 
to ^^^Au-t-^^^Au. In our earlier communication [sj, a 
prediction of E^ai for 238u^_238u ^g^g g^jg^ reported. It 
is worth mentioning that this is the first ever attempt 
that deals with as many as fifteen systems. Earlier at- 
tempts of mass dependence jS, 12, 16, 17, 26J had taken 
about six cases into account. Apart from the experimen- 
tal data, we also show our results for cr=35 and 40 mb. 
All curves are a fit of the form c.A^ using minimiza- 
tion procedure. The experimental data can be fitted by 
Texpt = —0.42079 ± 0.04594 whereas our present results 
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with cr=40 mb has T40 = -0.41540 ± 0.08166. The re- 
sults with cr=35 mb yields T35 = -0.43037 ± 0.08558. 
Exclusively, one can extract the following: 

• The present value of Texpt differs from the earlier 
reported results {Texpt. — —1/3 for A<200 5] and 
Texpt = -0.46 ± 0.06 _16| for A<400 mass). Note 
that in earlier attempts ^12, 16], the data of NSCL 
alone was used to fit the power law. In the present 
case, we have taken data from different sources 
amounting to 15 reactions. Therefore, the present 
analysis is universal in nature. 

• The present theoretical value r4o — —0.41540 ± 
0.08166 is the closest one obtained so far. In earlier 
reports 0, the Texpt was —1/3 whereas BUU model 
yielded -0.28<rth<-0.32. In another study [ll|, 
the Texpt was — 0.46±0.06 whereas BUU model has 
Tth — —0.41 ± 0.03. In other words, the present 
QMD model with a stiff EOS along with cr=35- 
40 mb can explain the data much better than any 
other theoretical calculations. Our present analysis 
has much wider mass spectrum than any early at- 
tempt. Some deviations in the middle order are also 
reported by other authors [^^ . The cr=40 explains 
E})ai in heavier nuclei whereas cr=35 reproduces a 
middle order nicely. 

• From the figure, it is also evident that a true cross- 
section for this energy domain should be between 35 
and 40 mb. This conclusion is very important since 
a wide range of masses was used for the present 
analysis. Our conclusion about the strength of the 
nn cross-section is in agreement with large number 
of earlier calculations on disappearance of flow and 
ot her p henomena in heayy-ion collisions |^, [lol l2p . 
milllilllillilslllllii. As noted in ref. jig, 
this value of nn cross-section is still much less than 
the actual averaged free nn cross-section which is 
about 60 mb 38]. In the case of fragmentation, a 
larg er cross-section of a—bb mb is also suggested 

HI. 

We have also tried to fit the balance energy in terms 
of other parameters such as the charge of the collid- 
ing nuclei. This attempt is shown in fig. 6 where 
Ebai is plotted as a function of the total atomic num- 
ber of the system. In the upper part, we display the 
full range of the systems whereas in the lower part, 
only heavier nuclei are taken into picture. A power law 
cxZ^ fits the data nicely. Now r is -0.46703 ± 0.04745 
(-0.45808 ± 0.08688) for experimental data (theoreti- 
cal results) which is larger compared to mass power law 
(r = -0.42079 ± 0.04594(-0.41540± 0.08166)) for ex- 
periment data (theoretical results). This difference in 
the slopes stems from the different charge/mass ratio in 
lighter and heavy nuclei. Interestingly, the value of r 
for heavier nuclei (see fig. 6(b)) is —0.59316 ± 0.0622 
(-0.58483 ±0.11409). This result, which is in agreement 



with ref. shows the dominance of the Coulomb inter- 
actions in heavier colliding nuclei. It would be of further 
interest to investigate whether the flow due to the colli- 
sion and mean field parts (at the balance energy) exhibit 
any mass dependence or not. We display, in fig. 7, the 
flow (pf"") at balance energy due to the collision (upper) 
and mean field parts (lower). Interestingly, we do not 
see any clear mass dependence. Rather very weak de- 
pendence (with T = -0.09427 ± 0.08379) exists on the 
system size. This is in agreement with ref. |32|, where a 
similar conclusion was drawn. 



IV. SUMMARY 

We have studied the mass dependence of the disap- 
pearance of flow in large number of colliding nuclei using 
QMD model. As many as fifteen reactions with masses 
between 47 and 476 were studied for the first time, where 
experimental balance energy was available. Our findings 
suggest a weak dependence of different cross-sections 
for heavier colliding nuclei in agreement with 12]. For 
the first time, fifteen different reactions were studied 
and our calculations with a stiff equation of state are 
in excellent agreement with experimental data. We 
could reproduce the slope of the power law [(xA^) very 
well. Our calculations suggest a cross-section of 35-40 
mb in this incident energy domain. We also showed 
that the collective flow due to mean field is attractive 
whereas it is repulsive for collision part. The balanc- 
ing of these two parts results in the disappearance of flow. 

This work is supported by the grant (No. SP/S2/K- 
21/96) from Department of Science and Technology, 
Government of India. 



6 



[1] W. Schcid, H. Miillcr, and W. Grcincr, Phys. Rev. Lett. 

32, 741 (1974). 
[2] B.A. Li and A.T. Sustich, Phys. Rev. Lett. 82, 5004 

(1999) . 

[3] Y.M. Zheng, CM. Ko, B.A. Li, and B. Zhang, Phys. Rev. 

Lett. 83, 2534 (1999). 
[4] J.J. MoUtoris and H. Stocker, Phys. Lett. B 162, 47 
(1985); V. de la Mota, F. Sebille, M. Farine, B. Re- 
maud, and P. Schuck, Phys. Rev. C 46, 677 (1992); G.F. 
Bertsch, W.G. Lynch, and M.B. Tsang, Phys. Lett. B 
189, 384 (1987); D. Krofcheck et ai, Phys. Rev. Lett. 
63, 2028 (1989). 
[5] G.D. Westfall et al, Phys. Rev. Lett. 71, 1986 (1993). 
[6] R. Pak, O. Bjarki, S.A. Hannuschke, R.A. Lacey, J. Lau- 
ret, W.J. Llope, A. Nadasen, N.T.B. Stone, A.M. Vander 
Molen, and G.D. Westfall, Phys. Rev. C 54, 2457 (1996); 
R. Pak et al, ibid. 53, R1469 (1996). 
[7] J.C. Angelique et al, Nucl. Phys. A614, 261 (1997). 
[8] J. P. Sullivan et al, Phys. Lett. B 249, 8 (1990). 
[9] D. Krofcheck et al, Phys. Rev. C 43, 350 (1991); C.A. 
Ogilvie et al., Phys. Rev. C 42, RIO (1990). 
[10] Z.Y. He et al., Nucl. Phys. A598, 248 (1996). 
[11] D. Cussol et al, Phys. Rev. C 65, 044604 (2002). 
[12] D.J. Magostro, W. Bauer, and G.D. Westfall, Phys. Rev. 

C 62, 041603(R) (2000). 
[13] G.D. Westfall, Nucl. Phys. A681, 343c (2001); R. Pak et 
al., Phys. Rev. Lett. 78, 1026 (1997); R. Pak et al., ibid. 
78, 1022 (1997). 
[14] D. Krofcheck et al, Phys. Rev. C 46, 1416 (1992). 
[15] W.M. Zhang et al, Phys. Rev. C 42, R491 (1990); M.D. 
Partlan et al, Phys. Rev. Lett. 75, 2100 (1995); P. Cro- 
chet et al, Nucl. Phys. A624, 755 (1997). 
[16] D.J. Magestro, W. Bauer, O. Bjarki, J.D. Crispin, M.L. 
Miller, M.B. Tonjes, A.M. Vander Molen, G.D. Westfall, 
R. Pak, and E. Norbeck, Phys. Rev. C 61, 021602(R) 

(2000) . 

[17] A. Buta et al, Nucl. Phys. A584, 397 (1995). 
[18] J. Singh, S. Kumar, and R.K. Puri, Phys. Rev. C 62, 
044617 (2000). 

[19] R.K. Puri and S. Kumar, Phys. Rev. C 57, 2744 (1998); 
Zhuxia Li, C. Hartnack, H. Stocker, and W. Greiner, ibid. 
44, 824 (1991). 

[20] C. Liewen, Z. Fengshou, and J. Genming, Phys. Rev. C 

58, 2283 (1998). 
[21] B.A. Li, Z. Ren, CM. Ko, and S.J. Yennello, Phys. Rev. 

Lett. 76, 4492 (1996). 
[22] H. Zhou, Z. Li, Y. Zhuo, and G. Mao, Nucl. Phys. A580, 

627 (1994). 

[23] D. Klakow, G. Welke, and W. Bauer, Phys. Rev. C 48, 
1982 (1993). 

[24] H.M. Xu, Phys. Rev. C 46, R389 (1992); H.M. Xu, Phys. 

Rev. Lett. 67, 2769 (1991). 
[25] B.A. Li, Phys. Rev. C 48, 2415 (1993). 
[26] H. Zhou, Z. Li, and Y. Zhuo, Phys. Rev. C 50, R2664 

(1994). 

[27] S. Soff, S.A. Bass, C. Hartnack, H. Stocker, and W. 

Greiner, Phys. Rev. C 51, 3320 (1995). 
[28] E. Lehmann, A. Faessler, J. Zipprich, R.K. Puri, and 

S.W. Huang, Z. Phys. A 355, 55 (1996). 
[29] S. Kumar, M.K. Sharma, R.K. Puri, K.P. Singh, and I.M. 

Govil, Phys. Rev. C 58, 3494 (1998). 



[30] A. Sood and R.K. Puri, Symposium on Nuclear Physics, 

45B, 288 (2002); A. Sood and R.K. Puri, VIII Int. Conf. 

on Nucleus-Nucleus Collisions, Moscow, Russia, June 17- 

21 (2003)- accepted. 
[31] A.D. Sood, R.K. Puri, and J. Aichelin, Phys. Rev. Lett., 

(2003)- submitted. 
[32] S. Kumar, R.K. Puri, and J. Aichelin, Phys. Rev. C 58, 

1618 (1998). 

[33] J. Aichelin, Phys. Rep. 202, 233 (1991); J. Aichelin and 
H. Stocker, Phys. Lett. B 176, 14 (1986). 

[34] C. Hartnack, R.K. Puri, J. Aichelin, J. Konopka, S.A. 
Bass, H. Stocker, and W. Greiner, Eur. Phys. J A 1, 151 

(1998). 

[35] H.W. Barz, J. P. Bondorf, D. Idier, and I.N. Mishustin, 

Phys. Lett. B 382, 343 (1996). 
[36] C. Roy et al, Z. Phys. A 358, 73 (1997). 
[37] B. Blattel, V. Koch, A. Lang, K. Weber, W. Gassing, and 

U. Mosel, Phys. Rev. C 43, 2728 (1991). 
[38] K. Chen, Z. Fraenkel, G. Friedlander, J.R. Grover, J.M. 

Miller, and Y. Shimamoto, Phys. Rev. 166, 949 (1968). 



7 



Figure Captions 

FIG. 1. The averaged {p^jA) as a function of 
^c.m.l^heam- Here we display the results at different 
incident energies using a stiff equation of state along 
with cr=40 mb. The reactions of 2°Ne+^"Al and 
4°Ar+45S(, j^j.g ]^5o fm/c whereas those of 64Zn+58Ni, 
93Nb+93Nb, i39La+i39La and is^Au+^^Uu are at 300 
fm/c. 

FIG. 2. The time evolution of (p*'') as a function of 
time. Here again results are for stiff equation of state 

along with cr=40 mb. 

FIG. 3. The (jJ*'") as a function of the incident 
energy. The results for different cross-sections of 55, 
40 and CMC are represented, respectively, by the solid 
circles, open squares and solid inverted triangles whereas 
for Cug and 20 mb are represented by solid diamonds 
and solid triangles. A stiff equation of state has been 
used. All lines are to guide the eyes. 

FIG. 4. The decomposition of (pf"") into collision 
and mean field parts as a function of incident beam 



energy. Here results arc displayed for (T=40 mb. Stars 
are the experimental balance energy. 

FIG. 5. The balance energy as a fimction of combined 
mass of the system. The experimental points along 
with error bars are displayed by solid stars whereas our 
theoretical calculations for (7=35 and 40 mb arc shown 
by open triangles and squares. The lines are the power 
law occ.A^. The solid, dashed and dash-dotted lines 
represent the power law fit for experimental points, with 
CT=40 and 35 mb, respectively. 

FIG. 6. (a) Balance energy as a hmction of atomic 
number Z. Here we display the experimental results 
along with our calculations for a=40 mb. The fits are 
obtained with minimization for power law function 
c.Z'^. (b) Same as (a), but for heavier nuclei. 

FIG. 7. The decomposition of the (p*'') at balance 
energy into collision and mean field parts. The results 
are obtained using a stiff equation of state along with 
(7=40 mb. 
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